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Abstract. The hyperensembles, introduced by Crooks in a context of non-equihbrium statistical 
physics, are considered here as a tool for systems in equilibrium. Simple examples like the ideal 
gas, the mean-field model, and the Ising interaction on small square lattices, are worked out to 
illustrate the concepts. 
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INTRODUCTION 



The starting point of superstatistics |lll,l2D is the assumption that the inverse temperature 
/3 of the Bohzmann-Gibbs distribution 

may itself be a fluctuating quantity, governed by a probability density /(/3). The result- 
ing probability distribution is then determined by 

/•oo 

Pf{x) = d^mqp{x). (2) 

The choice of the hyperdistribution /(/3) (called the entropic distribution in f^) depends 
on the application at hand. Recently, Crooks [3] suggested that this distribution should 
be determined by means of the Maximum Entropy Principle (MEP) because, in general, 
information is lacking about what choice of /(/3) is appropriate. Additional constraints, 
introduced when applying the MEP, then lead to a parametrised family of probability 
distributions. The latter he called a hyperensemble. 

The choice of constraints is crucial. Crooks proposes to use three constraints: the 
normalisation, the mean energy and the mean entropy. The result of applying the MEP 
is that the hyperdistribution is of the form 

/(/3) = 0j-f^P -e{H)p) . (3) 



Here, is the Lagrange multiplier controlling the mean energy, A controls the mean 
entropy. The canonical entropy is denoted S(j8) to distinguish it from the entropy of the 



hyperensemble. The normalisation is given by 



f'OO 

Z(0,A) = c(/3)d/3exp(A5(/3)-0(//)^). (4) 

The intention of Crooks is to use the notion of hyperensembles to study systems out 
of equilibrium. The present paper shows the usefulness of hyperensembles as a unifying 
concept for systems in equilibrium. In this context the parameter A plays no important 
role and may be taken equal to 1. 



THE IDEAL GAS 

Consider only the momenta pi,p2, - ■ ■ ,pn of an ideal gas containing A'^ particles. The 
Hamiltonian is 



(5) 



The canonical probability distribution is 



The entropy and energy equal 



5(/3) = -^ln/3, {H)p = '-^P-'- (7) 

After optimisation, the hyperdistribution becomes 

/(.) = i|e.p(-f,„,-e|). 

This is the inverse Gamma distribution, also called the inverse ;f^-distribution. It has a 
maximum at /3 = 0. Hence, the canonical probability distribution qp{p), with /3 = 0, is 
the most likely distribution at inverse hypertemperature 0. 

The ideal gas, together with the above hyperdistribution, was already considered in 
the context of superstatistics 10, 111]. 



THE MEAN-FIELD APPROXIMATION 

Consider spin variables ai, (727 ■ • ■ , (^n, taking on values ±1. The equilibrium states of 
the ideal paramagnet are the product states 

N 

^g(a) = ]^e3(i+^«)(l-e)5(i-^«), (9) 

n=l 
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FIGURE 1. Hyperdistribution fp_h{e) of the mean-field model with / = 1, = 5, /; = 0, and j3 = 0.6 
and 15 — 1.2, respectively. 



with < e < 1. The total magnetisation is (M)e = N{2£ — 1). The entropy is 

5(e) = -A^elne-A^(l-e)ln(l-e). (10) 
Here, we use £ as a parameter instead of /3 . Introduce the ferromagnetic Hamiltonian 

J N N 
^ = ~2 ^ Jm,n<ym(y„ - h Om. (11) 

m,n=l m=l 

Its average, using the equilibrium distribution of the paramagnet and assuming J^m = 0, 
equals {H)e = -\NJ{2e -if- hN{2e - 1), with 

1 ^ 

J ^ Jm,n- (12) 

m,n=l 

The optimised hyperdistribution with parameters /3 (instead of 6) and h is 



elne - (1 - e) ln(l - e) + ^/37(2£ - 1)2 + /3/z(2e - r 



(13) 



For /3 > /3mf = and h = the hyperdistribution has two maxima, corresponding 
to two equally likely states. See the Figure [T] At high temperature the paramagnetic 
state qeio) with £ = 1/2 is the most likely state. At low temperature the two states 
with e 7^ 1/2, maximising /j3^/j(e), are the most likely ones. This bifurcation from a 
paramagnetic state to a pair of ferromagnetic states as a function of the temperature ^ 
involves a spontaneous magnetisation of the system. 



In the quantum case the ideal paramagnet is described by a density matrix of the 
product form Pq'^ , with 

1 3 

Pe = ^ + Y.^kOk- (14) 
^ k=\ 

The O]^ are the Pauli matrices. The Hamiltonian of the Heisenberg model reads 

^ N 3 N 

^^~2 ^ Jm,nY,^'nk<^nk-hY,<^n3- (15) 

m,n=l k=l n=l 

The matrix is a copy of at site m. Without restriction, assume = 0. Then the 
average energy in the product state is (H) q = —NJ\ p +Nh63, with / as in the classical 
case. The von Neumann entropy of the product state equals 

Sip'f) = -f(l + |0|)ln^(l + |0|)-f(l-|0|)ln^(l-|0|). (16) 

One can now calculate the hyperdistribution fp,h{^)- Note that this is a classical distri- 
bution, not an operator. A maximum is attained when the following three equations are 
satisfied 

^|^l"T3^ = -2j^-^^'^-^^^'3, ^=1,2,3. (17) 

Assuming h ^ 0, the solution requires O2 = 63 = 0. The third parameter 63 must be a 
solution of the usual mean-field equation O3 — tanh(/3y03 + 63/1). 

THE MICROCANONICAL ENSEMBLE 

Following Boltzmann, the probability distribution of the microcanonical ensemble 
equals 

qE{x) = ^^d{E-H{x)), (18) 

with the density of states p{E) given by p(£') = j dxc{x)d{E — H{x)). The Boltzmann 
entropy is S{E) = fcfilnp(£'). The optimised hyperdistribution becomes 

/,W = ^exp(^W-e£)=|||).-«-. (19) 

The probability distribution pg (x) of the hyperensemble coincides with the canonical 
Boltzmann-Gibbs distribution 

poo 

pe{x) = / dEfe{E)qE{x) 
Jo 
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FIGURE 2. Number of configurations of a 2x2 and a 3x3 Ising lattice as a function of energy and 
magnetisation. On the vertical axis is the number of up spins m, on the horizontal axis the number of 
non-matching neighbour pairs «. 



Jo ^^m W) ^ ^^ ^^ 



= ^.-^H^). (20) 
Z(0) ^ 

The parameter 9 is the inverse temperature /3. 

As an example, involving a discrete phase space, let us consider the 2-dimensional 
Ising model on a square lattice of n x n Ising spins. Count the number of configurations 
C(m, n) with m spins up and n pairs of unequal nearest neighbours. See the Figure |2l 
The entropy is y{m,n) = lnC(m,n). The hyperdistribution is 

fp,h {m, n) = exp (2j87(A^ - n) - /3/i(A^ - 2m) ) . (21) 

At constant n and h = it is proportional to C{m,n). See the Figure [3l It shows two 
maxima, corresponding with the spin up and spin down phases. 



DISCUSSION 

The aim of the paper is to show that the concept of hyperensembles, as introduced by 
Crooks [3] in a context of non-equilibrium systems, fits well into the standard theory 
of equilibrium statistical physics. Several situations have been considered. In the ex- 
ample of the ideal gas the hyperensemble approach coincides with the superstatistical 
treatment of [3,[5l]. In the mean-field theory one starts from a simple model, such as the 




FIGURE 3. Normalised C{m,n) as a function of m at constant n = 16 for a 5x5 Ising lattice. 



ideal paramagnet, to solve more complex models under the constraint that the state of 
the system is an equilibrium state of the simple model. The present formulation in terms 
of hyperensembles is a mere reformulation of what is known since long. However, it can 
be applied in a very general context. Such an application beyond the traditional scope 
of mean-field theory can be found in [6], where a random walk model is used to model 
polymer behaviour. In the microcanonical ensemble the probability distribution of the 
hyperensemble coincides with the canonical Boltzmann-Gibbs distribution. The maxi- 
mum of the hyperdistribution may be degenerate. This is for instance the case in finite 
spin lattices with ferromagnetic Ising interaction. Hence, the most-likely microcanon- 
ical state is non-unique. Such a feature also occurs in mean-filed models and may be 
interpreted as a precursor of the phase transition occurring in the thermodynamic limit. 

REFERENCES 

1. C. Beck and E.G.D. Cohen, Superstatistics, Physica A 322, 267 - 275 (2003). 

2. C. Beck, Superstatistics: Theoretical concepts and physical applications, to appear in Anomalous 
transport: Foundations and Applications, ed. R. Klages et al. (Wiley, 2007), arXiv:0705.3832 

3. G.E. Crooks, Beyond Boltzmann-Gibbs statistics: Maximum entropy hyperensembles out-of-equi- 
librium, Phys. Rev. E 75, 041119 (2007), arXiv:cond-mat/0603120 

4. H. Touchette, Temperature fluctuations and mixtures of equilibrium states in the canonical ensemble, 
in Nonextensive Entropy, ed. M. Gell-Mann and C. Tsallis (Oxford University Press, 2004), pp. 159 - 
176. 

5. H. Touchette and C. Beck, Asymptotics of superstatistics, Phys. Rev. E 71, 016131 (2005). 

6. E. Van der Straeten and J. Naudts, The globule-coil transition in a mean fleld approach, 
arXiv : cond-m at/06 1 2256 



